Spin resolved Hall effect driven by spin-orbit coupling 
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Spin and electric Hall currents are calculated numerically in a two-dimensional mesoscopic system 
with Rashba and Dresselhaus spin-orbit coupling by means of the Landauer-Biittiker formalism. It 
is found that both electric and spin Hall currents circulate when two spin-orbit couplings coexist, 
while the electric Hall conductance vanishes if either one is absent. The electric and spin Hall 
conductances are suppressed in strong disorder, but survive in weak disorder. Physically it can be 
understood that the spinomotive transverse "force" generated by spin-orbit coupling is responsible 
for the formation of the spin Hall current and the lack of transverse reflection symmetry is the origin 
of the electric Hall current. 
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When a metallic sample is subjected to a perpendicu- 
lar external magnetic field, the Lorentz force acting on 
the charge carriers gives rise to a transverse voltage be- 
tween two edges of the sample, this is well known as 
the ordinary Hall effect. Actually the Hall effect fam- 
ily has numbers of important members. The anomalous 
Hall effect may occur even in the absence of an external 
magnetic field in a ferromagnetic metal with spin-orbit 
couplingi 2 i2i2i5i£iL& In the past few years it has been rec- 
ognized that the spin-orbit coupling may provide an effi- 
cient way to manipulate a spin resolved current in metals 
and semiconductorsi2ii2iiiii2ii&i£ii£ In a two dimensional 
electron gas (2DEG) lacking bulk and structure inversion 
symmetries, the effective Hamiltonian is given by 



H = 



P 

2m* 



+ X(o- x p y - a v p x ) + /3(o- x p x 



c V Vy) 



(1) 



where the second term is the Rashba spin-orbit coupling 
and the third one is the Dresselhaus spin-orbit coupling. 

(/i = x, y, z) are the Pauli matrices and the coupling 
parameters A and /3 have the dimension of velocity. Using 
the Heisenberg equation of motion the second derivative 
of the position operator r gives 



dt 2 



2m*(X 2 - (3 2 )c 



-p x z. 



(2) 



Compared with the Lorentz force brought by the mag- 
netic field upon a charged particle, the spin-orbit cou- 
pling produces a spinomotive transverse "force" on a 
moving electron. It has no classical counterpart as the 
coefficient is divided by h, but it reflects the tendency 
of spin asymmetric scattering of a moving electron sub- 
ject to the spin-orbit coupling. When charge carriers are 
driven by an external electric field, this "force" tends to 
form a transverse spin current. 

In this paper we calculate the spin and electric Hall 
conductances in a 2DEG mesoscopic system with Rashba 
and Dresselhaus coupling by using the Landauer-Biittiker 
formula and the Green's function technique. It is found 
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FIG. 1: Cross-shape device with four semi-infinite metallic 
leads. The spin-orbit coupling is supposed to exist in the 
shadowed area only, and the effect of the semi-infinite leads 
is treated exactly through self-energy terms. 



that both electric and spin Hall currents circulate while 
these two types of spin-orbit coupling coexist, but the 
electric Hall current disappears when either one is absent. 
The spin and electric Hall conductances are suppressed 
in strong disorder, but survive in weak disorder. The nu- 
merical results are in good agreement with the symmetry 
analysis of the system. 

We consider a cross-shape device with four semi-infinite 
metallic leads as shown in Fig. 1. The scattering region 
(shadowed part in Fig. 1) is described by the effective 
Hamiltonian in Eq. Q , and when it is treated as an L x L 
lattice with the tight-binding approximation, the model 
Hamiltonian reads, 

H = —t C i,cr C 3,0 + tso {^i,y ~ Jf,x) 

i 

where and tj? Q are the dimensionless parameters for 
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Rashba and Dresselhaus coupling strength in the unit of 
t, respectively, and the local spin current operator J? is 
defined aa~ 

J L = (4x CT ^ +Q ^ - h. c ) . (4) 

a, a' 

where a stands for the unit vector along axes of the lattice 
and /i stands for the direction of spin polarization^ 

The calculation of electric and spin currents is based 
on the Landauer-Buttiker formalisniiiSiiS Assume T"'g to 
be the spin-resolved transmission probability of electrons 
transmitted from spin channel (i of lead p to spin channel 
v of lead q, and V p to be the electric voltage in lead p, 
then respectively the outgoing electric current and spin 
current polarized along fi direction in lead p are 

ic p = (5) 

The transmission probability coefficients can be cal- 
culated by using the Green's function technique, 
T£g =Tr[r»G R T£G A ] . And the retarded and ad- 
vanced Green functions are given by G RA (E) = 
l/(E — H c — J2p=i^p' A )> where E is the electron en- 
ergy and H c is the model Hamiltonian in the shad- 
owed region in Fig.l. The retarded and advanced 
self energy terms introduced by the semi-infinite lead 
V, (n^W, Pj(T ' = -tJ2 m Xm,(Pi)e lkma Xrn{Pj)5 aa > and 
lip = (H R y where Xm(Pi) is the transverse mode wave 
function at site pi in lead p connected to the scat- 
tering region. It should be noted that in Eq.© fi 
may stand for an arbitrary direction of spin polariza- 
tion, and this is done by incorporating a transforma- 
tion in the definition of r p , that is, T£(j>i, a, pj, a') — 

ZtR^R'^Y^mXmipi) sin(k m a)x„i{Pj) and R is the ro- 
tation matrix to transform the eigenstates of a z to those 
of t • <j (f is a unit vector) iSl The Landauer-Buttiker for- 
malism has been applied extensively to study the spin 
transport in mesoscopic systems numerically^i*22i2i24 

In this paper we consider an initial electric current 
driven through leads 1 and 2, I\ = — /f, by setting the 
bias voltage V x = -V 2 = V /2 and V 3 = V A = 0. The cur- 
rents in leads 3 and 4 are perpendicular to the current 
through leads 1 and 2. Thus the electric and spin Hall 
conductances are defined as 

G cH = Ig/iVi-Vt); (7) 
Gsh = Ig/(Vi-V 2 ), (8) 

respectively, where the spin current has three compo- 
nents, p — x,y,z. Electric and spin Hall conductances 
are evaluated as functions of the Fermi energy for differ- 
ent ratios of Rashba and Dresselhaus coupling strength 
in Fig. [21 Generally speaking, the electric Hall con- 
ductance is symmetric about the Fermi energy Ef while 
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FIG. 2: Electric and spin Hall conductances as functions of 
the electron Fermi energy for different ratios of Rashba and 
Dresselhaus coupling strength. L = 40 for all the results in 
this figure. 



the spin Hall conductance is antisymmetric such that 
the spin Hall conductance vanishes at the band center, 
Ef = 0. This is consistent with the symmetry analysis 
for the tight binding Hamiltonian. 23 In the case of the 
pure Rashba or Dresselhaus coupling, the electric Hall 
conductance disappears, but the spin Hall conductance 
still exists. In the two cases of tf G = 1 and t§ Q = 1/2 
and of t R — 1/2 and t^ = 1, the electric Hall conduc- 
tances are equal. However the spin Hall conductances 
G z sH differ by a minus sign, with G x sH and G y sH swapped, 
and the former is consistent with Shen and Sinitsyn et 
al's works for free 2DEG systemsi££*2& A special case is 
at the symmetric point of tg Q = t R Q . The spin Hall con- 
ductance G z sH = 0, while G X H and G y sH are equal and 
non-zero, which means the current is polarized within 
the x-y plane. In this case the operator a x + a y com- 
mutes with the total Hamiltonian, and actually there 
is no spin flip in the scattering region^ On the other 
hand the longitudinal electric and spin conductances are 
also non-zero. The longitudinal conductances are about 
one order larger than the Hall conductances in magni- 
tude, i.e., I3/I2 ~ 10- And the electric conductance is 
also symmetric with respect to the Fermi energy, just 
like the electric Hall conductance, while the longitudinal 
spin current is antisymmetric. According to the symme- 
try properties of such a system^ under the transforma- 
tion: a x ^> a y , a y ^ a x , and a z — > — a z , tg Q — > t R Q and 
^so -» *SO> correspondingly G X H -> G v sII , G v sH -> G x sH , 
and G z sH —G Z H while the electric Hall conductance 
remains unchanged. At the symmetric point, t R Q = tg Q , 
it yields that G X H — G y sH and G z sH = 0. Our numerical 
results obviously agree with this symmetry analysis. 

The Hall conductances as functions of the Rashba cou- 
pling strength are also evaluated, with specific Dressel- 
haus coupling strength t® — 1.0 and electron Fermi en- 
ergy Ef = — 2.0t as shown in Fig. 03 It indicates clearly 
that the electric Hall conductance increases with increas- 







FIG. 3: Electric and spin Hall conductances as functions of 
the Rashba coupling strength with a fixed Dresselhaus cou- 
pling strength t®„ = 1.0 at E f = — 2.0t and L = 40. Similar 
results are obtained for the Hall conductances as functions of 
the Dresselhaus coupling strength with a fixed Rashba cou- 
pling strength, which are consistent with symmetry analyses. 



FIG. 4: Electric and spin Hall conductances divided by sam- 
ple size L as functions of the electron Fermi energy. Here the 
Rashaba and Dresselhaus coupling strength are equal and the 
results for L = 20, 30, 40, and 50 are shown simultaneously. 
It implies a size effect that both Hall conductances are pro- 
portional to L in this calculation. 



ing the Rashba couple strength, and reaches its maxi- 
mal value at tf Q ~ tf Q . Then it turns to decrease when 
*so > *soj an d approaches to zero for a large Rashba 
coupling strength. The figure shows that G z sH = and 
G X H = G v sH at tf a = tf as expected by the symmetry 
analysis. For a large spin-orbit coupling both electric 
and spin Hall conductance approaches to zero because 
the spin-orbit coupling in the scattering region forms a 
large potential barrier and the incident electrons will be 
completely reflected. Unlike bulk systems2^& where the 
spin Hall conductance in the clean limit has a universal 
value ±e/87r and the sign is given by the relative ratio of 
two coupling strength in Eq.JIJ, Fig. shows that the 
value of spin Hall conductance varies with the relative 
ratio as well as the sign, but the change of sign is com- 
patible with the bulk systems case. And this result is 
also compatible with the previous numerical work in the 
case of pure Rashba coupling^i2ii2& 

To see the finite size effect we calculate the electric 
and spin Hall conductances for L = 20, 30, 40, and 50. 
G c h/L and G s h/L as functions of Ef are plotted in Fig. 
21 We notice that these curves for different sizes fit a 
single one very well. Thus we conclude that both electric 
and spin Hall conductance are proportional to the size 
L of the sample. In other words, in our calculation the 
electric and spin Hall currents are determined by both the 
number of the incident channels and that of the outgoing 
channels. Thus the Hall currents induced by a specified 
longitudinal electric field are not proportional to the size 
L linearly, but to L x L. 

The disorder effect is an interesting issue in the spin 
Hall effect in 2DEG. It is still greatly controversial 
whether the spin Hall effect may survive when the im- 
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FIG. 5: Electric and spin Hall conductances as functions of 
logarithm of the disorder strength W/t. Results are obtained 
with equal Rashba and Dresselhaus coupling strength at Ef = 
—2.0t and L = 30. Standard deviations in the calculation are 
shown through the error bars. 

purity scattering is taken into account ,SL22i2i We con- 
sider the disorder effect by including the disorder poten- 
tial term Vdisorder = Z)i.<r=T,i e i c l.cr c i,° in Eq.© where 
e, are randomly distributed between [-W/2, +W/2]. Se- 
lectively the electric and spin Hall conductances, G c h 
and G x sH , for two couplings with equal strength are plot- 
ted in Fig. |3J G z sI1 is exactly equal to zero according 
to the symmetry. It shows that both electric and spin 
Hall conductances can survive in weak disorder, but be 
suppressed in strong disorder. We also did calculation 
for several other cases, and obtained similar results. The 
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case of pure Dresselhaus coupling is in agreement with 
Shcng et al's work for pure Rashba coupling. 23 

Physically the spin Hall conductance can be well un- 
derstood from the spinomotive transverse " force" caused 
by the spin-orbit coupling in Eq.©. The electric field 
drives electrons moving along the field such that the elec- 
trons with spin-up or -down experience opposite trans- 
verse "force" and thus a non-zero spin current is induced 
perpendicular to the field. The relative ratio of the two 
coupling strength determines the direction of the spin 
Hall current as the spinomotive force changes its sign 
around A = (3 and vanishes at the point. All calculated 
results are consistent with this. However, the spinomo- 
tive force is not a direct origin of the non-zero G c h, since 
G c h arises only when two couplings are present simulta- 
neously. From the symmetry properties of the system 
we notice that the Hamiltonian with pure Rashba cou- 
pling is invariant under the transformation: k x — » k x , 
k y — ► — k y and cr x — > —cr x , a v — > u v , a z — > —a z . 
Similarly the Hamiltonian with pure Dresselhaus cou- 
pling is invariant under the transformation: k x — > k x , 
k y — ► — k y and a x — > cr x , a v — > —a y , a z — > — a z . This is 
why the electric Hall current vanishes in these two cases, 
while the spin Hall current circulates because there is 
no symmetry constraint on it as both k y and a z change 
their signs under such transformation. On the other 
hand, the Hamiltonian with both Rashba and Dressel- 
haus couplings does not possess the reflection symme- 



try of k x — > k x , k y — > — k y . Therefore the coexistence 
of both couplings breaks the reflection symmetry of the 
system, which makes the electric current not parallel to 
the electric field such that it gives rise to a nonvanishing 
Hall conductance G c h- This unconventional Hall con- 
ductance may be related to some discussions in terms 
of the anomalous Hall effect due to parity anomaly and 
additional band crossings. Moreover, since the diagonal 
spin conductance is non-zero in this case2&, the diagonal 
spin current along leads 1 and 2 might generate a charge 
Hall current via the reciprocal spin Hall effect ^ 

In conclusion, we studied the electric Hall conductance 
as well as the spin Hall conductance for a finite-size sys- 
tem with four leads. Both electric and spin Hall con- 
ductances are non-zero when both Rashba and Dressel- 
haus coupling are present, thus the current is actually 
spin polarized. Unlike the anomalous Hall effect, the 
present electric Hall current is driven by the spin-orbit 
coupling, not by the exchange coupling with the magnetic 
impurities^ This effect also differs from the one resulted 
from a spin polarized current via the Rashba coupling. 34 
Though the incident current is not spin polarized, the 
Hall current is polarized in our case. 
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